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Abstract 

In a small window of phase space, chiral perturbation theory can be used to make 
standard model predictions for tau decays into two and three pions. For r — > 
27ri/T-, we give the analytical result for the relevant form factor Fy up to two loops, 
then calculate the differential spectrum and compare with available data. For r 
SirUr, we have calculated the hadronic matrix element to one loop. We discuss the 
decomposition of the three pion states into partition states and we give detailed 
predictions for the decay in terms of structure functions. We also compare with low 
energy predictions of meson dominance models. Overall, we find good agreement, 
but also some interesting discrepancies, which might have consequences beyond the 
limit of validity of chiral perturbation theory. 



1 Introduction 



Semileptonic decays of the heavy tau lepton into a tau neutrino and a hadronic 
system offer a unique laboratory to study the standard model and especially low- 
energy QCD. A variety of multi- meson final states with invariant masses from the 
production threshold up to the tau mass of about 1.8 GeV can be studied. These final 
states can also be produced using hadronic initial states (such as pion-nucleon or 
nucleon-nucleon collisions). The production in tau decays, however, is advantageous 
in that the initial state is simple, clean and well understood. Whereas some of the 
final states (e.g. two or four pions) can also be produced using an electromagnetic 
current, i.e. electron-positron annihilation, other states (e.g. three pions in an / = 1 
state) can only be produced through the weak current. And in the case of the states 
which can be produced both in tau decays and in electron-positron annihilation, 
such as the two and four pion final states, tau days now compete well in statistics 
with electromagnetic production. 

For inclusive semileptonic tau decays, the tau mass of about 1.8 GeV is large 
enough to allow the application of perturbative QCD and in fact offers a unique 
possibility to measure the strong coupling constant as{^) at the low scale fi = M^. 
|l[]. In the case of exclusive semileptonic decays, calculations based on a systematic 
use of the QCD Lagrangian are not available up to now. In these decays, the probe 
testing the hadronic current carries a momentum Q which is mainly in the interme- 
diate energy region, the most difficult one for the study of strong interactions. In 
this region theoretical predictions have been obtained by using some kind of phe- 
nomenological models or approximate methods, such as quark models p|, ^, vector 
meson dominance ^, ^, |^, ^ |^, tree level calculations from effective Lagrangians 
10[| , or unitarization of current algebra results |]lT] . 

A small fraction of these decays, however, happens with very low Q^, below 
the mass of the lightest resonance. In this region the only active fields are the 
pseudoscalar mesons and one can use an effective Lagrangian to describe their in- 
teractions. This effective theory, called Chiral Perturbation Theory (CHPT), is a 
systematic method to calculate QCD matrix elements at low energy by means of an 
expansion in powers of the external hadronic momenta |jT2|, |13|, |14|. 

The limits of applicability of CHPT do not allow to give predictions for integrated 
decay rates, which would involve up to M^. Furthermore, tau decays into more 
than a single pion are dominated by resonant intermediate states, such as the p(770) 
in the two pion channel and the ai(1260) in the three pion channel. And so there 
are only relatively few events with small Q^. For these reasons, in the past, CHPT 
has been considered as not very interesting for tau decays. 

However, with the present high luminosity machines such as LEP and CESR, 
and even more so with future b and perhaps tau-charm factories, tau physics has 



2 



turned into an era of precision measurements, exploring very small branching ratios 
and studying details of differential distributions. Thus the small regime which 
is interesting for CHPT is now becoming accessible. 

We would like to mention that there is another small corner of phase space 
where a different systematic expansion of the hadronic current becomes possible. 
Heavy meson chiral perturbation theory [|15| can be applied to tau decays into a 
vector meson and a pion (such as r — >■ pyrz/,-, r — >■ K*^!/^-), if the momentum of 



the pion is small in the vector meson rest frame |16]. A complete calculation which 



includes vector meson decay and interference effects between different vector meson 
amplitudes, however, is still missing in this approach. 

Another reason why CHPT is relevant to tau decays is the fact that it can be 
used to test phenomenological models or fix some of their parameters. Indeed, the 
0{p^) prediction of CHPT in the limit of vanishing quark masses has been used to 
normalize vector meson dominance models in ^ |^, R, Q . In the present paper we 
will extend the CHPT prediction to higher order in p , and it will be a severe test 
for models if they correctly reproduce these higher orders. 

The expansion parameter of CHPT is (5^/(47rF^)^, with AtiFt^ = 1.2 GeV, so we 

are interested in \fQ^ below 500 . . . 600 MeV. Hadronic final states with a single 
pion or kaon can be predicted directly from and and so there is nothing 

interesting CHPT could teach us here. Final states with two and three pions allow 
for a reasonably large region of between threshold and the limit of applicability of 
CHPT. Already with four pions this region has almost disappeared. Moreover, the 
phase space for a n pion hadronic state, with close to threshold — > {nM^^Y, 
opens proportional to 

r— N (3n-5)/2 



(see Sec. 3.1 below). The exponent is 1/2 for two pions, 2 for three pions and 7/2 for 
four pions. So in the case of four pions the small interesting region for CHPT is even 
more suppressed by the phase space. As for final states with kaons, the threshold 
for a pion-kaon state is + Mk = 634 MeV . Furthermore the i^*(892) resonance 
is very close 



Therefore we find reasonable to try a CHPT calculation only for the 2tt and 3tt 
final states and these are the ones we will discuss in this paper. The two pion state 
is determined by only one form factor, namely the vector pion form factor. This 
can be measured in a few other processes, like e~^e~ —>■ 27t or vre scattering. Tau 
decays can provide an interesting cross check measurement, and we will investigate 
whether tau decays may become competitive in statistics. The three pion state, 
however, can only be produced in tau decay, and has a very rich and interesting 
substructure, which we will study in detail. 

Our paper is organized as follows: A brief review of chiral perturbation theory 
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is given in Sec. 2. In Sec. 3, we discuss the general structure of the phase space, 
the hadronic matrix element for two and three pions, and the differential decay 
rate. In this section we also discuss general properties of the three pion final state, 
regarding isospin invariance, the classification in terms of partition states and the 
definition of structure functions. We then calculate the hadronic matrix elements in 
CHPT in Sec. 4. Sec. 5 is dedicated to our numerical results, and in Sec. 6 we state 
our conclusions. In the Appendix we summarize our main conventions, display the 
decomposition of the three pion form factors into partition states and collect the 
main formulae regarding the definition of the structure functions. 



2 Chiral Perturbation Theory 

The relevant matrix elements for r decays into pions are of the form: 

(7r^^(pi)---^^"(p>ut|/;(0)|0) (2) 

1 1 

where = Vj^ = ^^T^^lt^Q when n is even, and -^^ = = ^'^T^lfil^Q when n 

is odd. A by now standard method to calculate such matrix elements in QCD at 
low energy is Chiral Perturbation Theory, CHPT. In this framework one uses an 
effective Lagrangian that respects the chiral symmetry properties of QCD, and that 
has the pions as the only "active" fields. Of course this Lagrangian is expected to 
be valid only up to energies which are well below the threshold for the production of 
heavier hadronic states. For more details about this method we refer the reader to 



the fundamental paper by Gasser and Leutwyler [|T3| and to a number of excellent 
reviews which are presently available in the literature |T^. Here we simply sketch 
the basic ideas and introduce the relevant notation. 

We consider the effective Lagrangian relative to two flavors in the isospin limit 
m = rriu = irtd. This Lagrangian contains an infinite number of terms; however it 
can be expanded in powers of derivatives and quark masses. One power of the quark 
mass will be counted as two powers of derivatives^ One will then have: 



^cff = £,2 + hCi + n'Ce + ... . (3) 

The leading order Lagrangian starts at 0{p^) and is the nonlinear a-model La- 
grangian in the presence of external fields, which we represent here in matrix form. 



^This means we are applying the "standard" CHPT counting rule. For a different counting rule, 
leading to a different ordering in the effective Lagrangian (the so-called "generalized" CHPT) , see 
Ref. Ml. 



4 




£2 = — (D^'t/D^f/^ + M^(f/ + t/t)) 

Ml = M^[l + 0{m)] 
= 2Bm 

= F[l + 0(m)] (4) 

B is proportional to the quark condensate (0|Mn|0) and the unitary 2x2 matrix U 
contains the pion fields, 

U = 

F F^ 

(5) 

The external fields v'', and af,, which we have introduced, have to be treated as 
external sources for the quark currents and A'^ respectively; in other words, 
the currents coupled to and in the effective Lagrangian are the low energy 
representation of the quark currents. In this framework these currents are expanded 
in powers of derivatives and quark masses, and are nonlinear in the pion fields, e.g. 
the axial vector current, 

K = '-^{r\U^D,U - UD,U^)) + 0{p^) = [-Fd,<p' + 0(0^)] + 0{p') . (6) 

Despite the fact that this Lagrangian is nonrenormalizable, one can use it to 
calculate matrix elements with the standard perturbation theory. As we have em- 
phasized in Eq. (^) the expansion parameter is h. This automatically produces an 
expansion of the matrix elements in powers of momenta and quark masses. As for 
the matrix elements in question, tree diagrams from £2 generate leading order con- 
tributions I, |, 1^, H, while one-loop diagrams yield terms at next-to-leading order. 
The occurring divergences in the loop contributions {in d = 4 dimensions) can be 
absorbed by introducing the effective Lagrangian at O(p^) [p!3| . 
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i 


k 


It 


1 


-1.7 ± 1.0 


1/3 


2 


6.1 ± 0.5 


2/3 


3 


2.9 ± 2.4 


-1/2 


4 


4.3 ± 0.9 


2 


6 


16.5 ± 1.1 


-1/3 



Table 1: The set of /j and corresponding 7j which we need for the calculation of the 
matrix elements in question. The values are taken from [EO] for i = 1,2, and from 



13] for all the others. The 7^ determine the relation between the k and the /[(/i). 



F^'^L = d^'{v''±a^)-d''{v^'±a^')-t[v>'±a'',v^'±a''] (7) 

where we omitted terms which contain external fields only. Since we disregard singlet 
vector and axial vector currents (i.e. (v^) = (a^) = 0) there is no contribution from 



the anomaly at 0(p ) ||T3 



The coupling constants k are split in a divergent and a finite piece and are 
scale-independent by definition, 

^"1-^^ = 7.^ + /i;f/[(/^) + O(rf-4) = (8) 
d/i Idtt^ ct/x 

The finite parts have been determined phenomenologically for the first time by 
Gasser and Leutwyler ]13[. We adopt their notation and use instead of the /[(/i), 
the finite and scale- independent quantities k, 

'' = (3^^) '"TIT W 

The up to date values for the relevant /j and the corresponding 7j are listed in Table 
|l|. For completeness we give the expressions for the pion decay constant and the 
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pion mass up to and including O(p^), 



Ml 



l + T7r^k + 0{M^) 



M' 



167r2F2 
M2 



327r2F2 



[3 + 0(M^ 



The mass splitting M^± — M^o is proportional to 



(10) 

m^Y and thus may be 
139.57 MeV andQ = 



neglected. In the numerical evaluation, we will use 
93.1 MeV. 

If one wants to go beyond the next-to-leading order, one has to calculate two 
loop diagrams with C2, and one loop diagrams with one vertex from £4. Again 
these diagrams will be divergent, but this is not a problem since at the same order 
one has contributions from tree diagrams with the Cq Lagrangian (which has been 
constructed in the case of three lights flavors l^ll). By defining appropriately the 
new coupling constants occurring in this Lagrangian, one is able to remove the 
divergences at the next-to-next-to-leading order, and get finite matrix elements. 
In order to have numerical predictions at this level one has to find a way to pin 
down or at least to estimate the finite parts of the new low energy constants. At 
the moment this has not been done yet in a systematic way: in Sec. [4.1| , when 
calculating the matrix element of the r decay into two pions to two loops, we will 
show how in a specific case one can try to circumvent this problem. 



3 Phase Space, Matrix Elements and Structure 
Functions 

3.1 Phase Space Considerations 

The two-pion phase space is given by 

$2. = ttt::^ / ^(1-^1(1-^1 fill 



1287r3 y M2 \ Ml \ g2 



(the reader is referred to App. |A.1| for our conventions) where the integral is from 
threshold 4M^ up to M^. Close to the two pion threshold, this implies 

1 1 /i_^)(J^_2M.)-'^ (12) 



dQ2 1287r3My' V ^^r 



^This means that we are neglecting 0(a) corrections to the decay width r(7r i^-v^^ from 
which Ftt is extracted, for more details see Ref. 1231. 
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The three pion phase space is 



~ 20487r5 j^Mi M2 1^ ) Jam^ Si 

X [{s, -Q'- Mlf - AQ^MlY'^{s,{si - AM^Y'^ (13) 
Close to threshold 9M^ this can be approximated by 



dQ'^ 2i033/27r4 I ]Vf2 



3M,j (14) 

By induction we can show that the phase space for a n pion hadronic state, 
with close to threshold, — > (riM^)^, opens proportional to: 

^«(v/^-nM.)'*-="^ (15) 

For n = 2, the exponent is 1/2, for n = 3 it is 2, recovering the above results. 

It is clear that in general, the more pions there are, the slower the phase space 
opens at threshold. Therefore for four and more pions, it is not only the high 
invariant hadronic mass, but also the behavior of the phase space at threshold which 
prevents the application of CHPT. 

3.2 Two Pion Differential Decay Rate 

The hadronic matrix element of the decay into two pions is characterized by a single 
form factor, Fy(Q^), 

H>^ = (7r-(pi)7r°(j92)|F^ - A>^\0) = V2{p, -p^YFy {Q^) (16) 

Only \Fy{Q'^) \ can be measured, and it can be obtained by measuring the differential 
distribution in using 

Here we have normalized to the electronic branching ratio of the tau: 



1927r3 
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3.3 Structure Functions and the Three Pion Differential De- 
cay Rate 

3.3.1 Form Factors and Isospin Relations 

The most general form of the hadronic matrix elements for the tau decays into the 

2-K~'K^ and 2ti^-k~ final states, compatible with the requirements of Lorentz, isospin 
and G parity invariance and Bose symmetry is given in terms of three functions 
F, G and H which have to satisfy the following property: 

F{S2,SI,S^) = +F(si,S2,S3) 
G(S2,S1,S3) = +^(51,52, S3) 

H{s2,Si,Sz) = -H{si,S2,S^) . (19) 

The matrix elements are then given by 
(vr°(pi)7r°(p2)7r-(p3)|^;(0)|0) 

= G'(si,S2,S3)(pi+P2),i + -H"(si,S2,S3)(pi -P2)m + -^(si,S2,S3)P3^ (20) 

(7r-(pi)7r-(p2)7r+(p3)|^;(0)|0) 
= G^'^Hsi>«2,S3)(Pi +P2)n + H^'^\si,S2,Sz){px - P2) n + F'-^\si, S2, s^)pz ^ . 
The form factors for 2ti~ti^ and 2ti^ti~ are related by isospin symmetry 

F^'^\si, S2, S3) = [^(^2, S3, Si) + G{S3, Si, S2) - H{S2, S3, Si) + H{S3, Si, S2)] 

G''^\si, S2, S3) = ^ [Fis2, S3, Si) + F{S3, Si, S2) + G{S2, S3, Si) + G(S3, Si, S2) 

+H{S2, S3, Si) - H{S3, Si, S2)] 

H^~^\si, S2, S3) = ^ [F(S2, S3, Si) - F(S3, Si, S2) - G'(S2, S3, Si) + G'(S3, Si, S2) 

-if(s2,S3,Si) - if(s3,Si,S2)] (21) 

Alternatively one can use a decomposition of the matrix element into only two 
functions, -Fi(si, S2, S3) and -^^(si, S2, S3), with Fs even under the exchange of the 
first two arguments, and Fi of mixed behavior: 

= [-Fl(si,S2,S3)(pi - PsT + Fi{s2, Si, S3){p2 - P3Y]Tf,u + FsQi^ 

9 



{n~ip,)n-{p2)n+ip3)\A-^i0m 
= S2, ss)ipi - p,r + f[-^\s2, Si, s,){p2 - p,r]T,, + F^^^Q, , (22) 

where 



Q fiQi/ 



Qt, = {P1+P2+P3),, 

= S1 + S2 + S3- 3M^ . (23) 

The decomposition into Fi and Fs has the advantage that these form factors corre- 
spond to a definite overall spin (viz. Fi corresponds to spin 1 and Fs to spin 0), and 
therefore the structure functions (see next Sec. ^^ ) are usually expressed through 
Fi and Fs- If F, G and H are known, we can calculate Fi, Fs through 

FsiQ'^, 81,82) = aF{si,S2,Ss) + (l- a)G{si, 82,83) - pH{si, 82,33) 
si + 82- 2Ml 



a 



2Q 



2 



A completely analogous relation holds for the form factors of the all charged matrix 
element. 

Let us emphasize two facts regarding the two charge modes. Firstly, the two 
matrix elements for the final states 2n^n^ and 2n'^n^ are not independent. If one 
knows the matrix element for one of the two states, the other one can be calculated 
using isospin symmetry. Secondly, however, isospin symmetry does not require that 
the form factors and decay rates are equal for the two modes. This fact can be seen 
if we decompose the three pion states in terms of partitions [p^] . 



3.3.2 Classification in Terms of Partitions 



In Ref. p4[. Pais introduced a classification of pion states with overall isospin 
/ = or 1 in terms of correlation quantum numbers [A'^iA^2A'3]. The three integer 
quantum numbers Ni are partitions of the total number of pions 



A^i > A'2 > A's > 
Ni + N2 + N3 = N 



(25) 
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Each state [A''iA'^2^3] is characterized by its symmetry property under the exchange 
of some of the momenta pi, ■ ■ ■ ,Pn- Such a state is easily constructed with the help 
of a Young tableau: each Young diagram must have three rows with A^^i, N2 and 

cells in the first, second and third row, respectively. The cells must then be 
filled with numbers going from 1 to A^, with the only rule that all the numbers in 
the rows (columns) must be organized in increasing order from left to right (top 
to bottom). The rule to construct a pion state from a tableau is very simple: one 
has to symmetrize with respect to the exchange of the momenta with the indices 
which are in a row, and antisymmetrize with respect to the momenta with the 
indices which are in a column. The order of these operations of symmetrization or 
antisymmetrization is not important, but must be fixed once and for all. 

Remarkably, all the states belonging to the same class defined by the partition 
[N1N2N3] share some common properties about isospin and charge distributions: 

1. the overall isospin / is uniquely determined and it is / = if A^^i — and 
N2 — N3 are both even, and 1 = 1 otherwise; 

2. the states in a class [A^iA^2^3] are composed by A^3 subsystems of three pions 
with 1 = and N2 — N3 subsystems of two pions with / = 1, and A''i — N2 
remaining single pions (trivially SN^ + 2(A^2 — ^3) + (^1 — ^2) = N); 

3. the N pion states which we are describing contain all possible charge distribu- 
tions (e.g. for N = 2 and zero total charge, they would contain both 7t~^tt~ and 

states). The probability of a state to contain a given charge distribution 
is a "class property" , i.e. is uniquely determined by the partition [N1N2N3] to 
which it belongs. 

For a more detailed account of the properties of these N pion states we refer the 



reader to the original article by Pais ||2J]. We now concentrate on the case of our 
interest = 3. 

In this case we have three possible partitions: [300], [210] and [111]. The [111] 
corresponds to 7r''"7r°7r~ in an overall 1 = state (e.g. from the decay to 37r). 
The remaining two partitions [210] and [300] have 1 = 1 and so they can occur in 
r SnUr- 

These two partitions differ in their branching ratios into the two charge distri- 
bution states. The [210] decays equally into 27r^7r^ and 27r°7r^: 



BR(27r 



TT 



whereas the [300] state prefers the all charged mode 



[210] : = 1 (26) 



, , BR(27r-7r+) , , 
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(^,^) immediately lead to the inequalities obtained in Ref. |^ 



1 ^ BR(27r07r-) ^1 1 ^ BR(27r-7r+) ^ 4 

5 - BR(all (Svr)-) " 2 ' 2 " BR(all (Svr)-) " 5 ' 

Experimentally, the branching ratios of the r into the two states are equal within 
the errors, so that certainly the [210] strongly dominates and a possible small ad- 
mixture of the [300] state (which, as we will show below, is required by CHPT), 
has not yet been established. Note that if the decay occurs only via a decay chain 
T ^ aii^T-, ai ^ p7T and p vrvr, as in vector meson dominance models, there is 
only the [210] state (because of the p resonance, there is one two-pion subsystem 
with 7 = 1, i.e. N2 — = 1), and both decay charge modes are produced with 
equal rates. 

If one has an analytic expression for the form factors in each of the two charge 
modes, one can easily construct the matrix element of a given partition state, by 
following the rules which we described above. The complete decomposition of each 
of the form factors F, G and H into the three partition states (two states for the 
partition [210] and one for [300]) is described in App. |A.2| . Here we give in an 
obvious notation the decomposition only of Fi, since it is the most important form 
factor in the numerical analysis: 



.[300] p[210] 



Fi = 

= 2Fr''^-Ff'"J , (2J 



where with we have indicated the sum of the two states belonging to the [210] 



class. From this decomposition the branching ratios given in Eqs. ( |26|j27| ) follow, 
and since F^^^^ = _p(+) ^ jg qIqq^y that the [300] partition is absent in the 
VMD model. 

3.3.3 Structure Functions and Differential Decay Rate 

The differential decay rate for a general hadronic decay is determined by 

dT{T STTUr) = cos^ Qc^y.vH^'' dP S^^^ (29) 

where the hadronic and leptonic tensors are 

V:=L^(L,)t, H^,:=H^{H,)\ (30) 

where = (hadronic final state|yl~(0)|0). The decay is most easily analyzed in 
the hadronic rest frame, and we can write 

L^,H>^' = Y,LxWx (31) 

X 
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In general, H^'^ can be characterized by 16 independent real functions. In our case 
of a three pion final states, there are restrictions due to G parity and Bose symme- 
try, which leave 9 independent structure functions Wx- These hadronic structure 
functions Wx depend on the kinematics only through the hadronic invariants Si, S2 
and Q^. The angular dependence is contained fully in the corresponding leptonic 
Lx- For details, see App. [A.3| below and p6|, 



There are four structure functions, Wa, Wc, Wd, and We, which arise from the 
spin-1 part of the hadronic current, i.e. they depend on Fi{Q'^ , si, s2). A single 
structure function, Wsa, arises from the spin-0 part and depends on Fs, and four 
functions, Wsb, Wsc, Wsd and Wse are due to interference between spin-1 and 
spin-0 amplitudes. 

The structure functions can be measured by observing angular distributions and 
taking moments (m) with respect to products of trigonometric functions of these 
angles. 

In the numerical evaluation in Sec. ^]2] we will plot Si, S2 integrated structure 
functions Wx 

WA,C,SA,SB,SciQ'^) = J dSidS2WA,C,SA,SB,SciQ'^,Sl,S2) 

Wd,E,SD,Se{Q'^) = J dSids2Sign{si- S2)WD,E,SD,SEiQ'^,Si,S2) (32) 

Without the energy ordering sign(si — S2), the relevant wx would vanish due to Bose 
symmetry. 

The integrated decay rate is determined by Wa and Wsa only, the other functions 
give vanishing contributions after integration over the angles. We have 

2. 1 dQ'Ml-Q^jl .^2^^1^^2Qn ) 

4 Calculation of the Hadronic Matrix Elements 
4.1 Two Pion Decay 

The hadronic matrix element which is relevant for the r decay into two pions is the 
following: 

(7r^(pi)7r'(p2)out|V;^(0)|0) = te'^^p, -p2),Fy{s) , (34) 

where s = Q'^ = {pi + P2Y, and = -^q^y^T^q- In the framework of CHPT, Fv{ s ) 
was calculated by Gasser and Leutwyler [O] to one loop, and by Gasser and Meifiner 



| 28[| up to two loops, by using a three times subtracted dispersive representation. The 
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function to be integrated inside the dispersive integral is a particular combination 
of the vector form factor and the / = 1, P-wave vrvr scattering amplitude at tree and 
one loop level. In Ref. PS[ the dispersive integral was calculated numerically. Here 
we are able to give a compact analytic expression of this integral 0: 



(35) 



Ml 



X 



167r2F2 [ 9 
/ Ml V 



- 1 + 247rVV(a;) - — 



X 



60 



IGn'^Fl 



+ ! + 
2 X 



1^(1+ 24vrVV-(x; 



30 



+ 
+ 



3191 



■ X 



223 
216 



16 vr^ X 



■ X 



9 



540 



37x + 15 



6480 

^ ( 7 x^' - 151 X + 99 ) J{x) + ^ ' 



(36) 



X 



30a;2 + 78 s- 128)7^1 (x) 
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+ 8n' [x'-—x-2j Ki{x) | . 

Where we have used the following functions: 
1 



J{x' 
Ki{x] 
KAx 



167r2 
1 F^ix 



F{x) + 2 



2 ' 



167r2 a 
1 Fix) 1 1 



167r2 xa"^ 327r2 xcr^ 



FHx) 



a' 



(37) 



+ 



with 



487r2 xcr^ Ixo"^ 



F(x: 



F^(x) + Tx^a^F{x 



TT 



192 327r2 ' 



a In ■ 



CT-1 

a+ 1 



(38) 



0" 



1 -4/a; 



The functions -^j(x) were recently introduced by Knecht et al. and like 

J( X ) are analytic everywhere apart from a branch cut from 4 to oo, and go to zero 
as X — i> 0. 



^The integrals one has to calculate here are similar to the ones that occur in the tttt scattering 
amplitude to two loops, see Ref. pol PQ| 
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All the constants which occur in /y ( s/M"^ ) are known (see Sec. 0). The subtrac- 
tion constants {r'^)v and Cy are calculable in CHPT and can be expressed in terms 
of the low energy constants /[(/u), chiral logs, and the new low energy constants 
which appear in Cq. With this representation of the subtraction constants given by 
CHPT one automatically satisfies the relevant Ward identities, up to the order at 
which one is working. We do not give this explicit representation here because up 
to now there is no information on the numerical value of the new low energy 
constants. In the future, with more accurate data on various low energy processes, 
and more two loop calculation available, one could try to pin down at least some 
of them, but this will require a considerable amount of work and it is beyond the 
scope of our analysis. 

We adopt in the following the notation of Gasser and MeiBner [EH] and write 



167r2F2 
1 

mm 



Ml 



+ 



167r2F2 
1 



/i 



60M2 167r2F2 



/2 



+ 0{Ml) 



(39) 



/i and /2 contain all the contributions at the two loop level. In |3T| the pion charge 
radius squared {r'^)y has been determined from experimental data by means of a 
simple model for Fy which contains (r^)y as the only free parameter. The authors 
obtain the result (r^)y = 0.431 ± 0.010 fm^, including the systematic error as a 
constraint in the normalization of the data. We repeat the fit with the expression 
(^) leaving (r^)y and Cy as free parameters. Furthermore we include a theoretical 
error, leading to 



0.431 ± 0.020 ± 0.016 fm^ 
3.2 ±0.5 ±0.9 GeV"^ 



(40) 



where the first and second errors indicate the statistical and theoretical uncertainties, 
respectively. We reproduce the central value of the radius squared given in |3l| with 
a larger statistical error, because we fit two parameters simultaneously: If we keep 

reduces by factor of two. The central value 



:y fixed, the statistical error in 



of is rather close to the value obtained by resonance saturation. 



4.1 GeV" 
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However we observe that for both parameters the theoretical uncertainties are 
of the same order of magnitude as the statistical errors. Unless one has a way to 
keep these theoretical uncertainties under control, we do not see how {r'^Yy can be 
determined with the accuracy indicated in |pT|. Note that we are not able to fix the 



low energy constant from our fit, since we do not have independent informations 
on /i. In the numerical evaluation, we will use the values in (^0]). 
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4.2 Three Pion Decays 

As we have seen in Sect. p.3.1| , it is sufficient to discuss only one of the two matrix 

elements, the other one can be calculated via the isospin relations (^T]) . So we will 

consider only the one with two neutral pions. 

First of all a general consideration: this matrix element contains a pole term in 
due to the direct coupling of the axial current to the pion. So the matrix element 

can be written in general as: 

+ H{si,S2,Ss){pi-p2)^, + F{si,S2,S3)p3f, , (41) 

where A^ttI-s, t, m) is the vrvr scattering amplitude as defined, e.g. in Ref. [1^. Note 
that the separation between the pole term and the barred form factors F and G 
is not unique, however one can split them such that the coefficient of the pole is 
exactly the tttt scattering amplitude, and therefore define in this way F and G. 

The calculation of the form factors is done by expanding them in powers of 
momenta and quark masses: 

R = ^^(^^°^ + |^ + ---) R = F,G,H 

A.. = + + , (42) 

where the superscript (n) indicates a contribution of order p". (We remark here that 
a tree diagram from the Lagrangian Cn gives a contribution of order to scattering 
amplitudes but of order to form factors.) 

At tree level if = for the simple reason that its antisymmetry under exchange 
of Si and S2 cannot be satisfied with a constant. As for F^^^ and G^^\ a constant 
satisfies their symmetry properties (0), and we find 

G^") = 1 . (43) 
At the one loop level we have the following results: 

F^'^ = lM'Mh) + Jm-^{s^-S2)[J{h)-J{h)]-ls3J{^^^ 
^ ' 2h{s3 - 2Ml) + [2(^1 + S2 + S3- AMI) - QhMl 



967r2 



Msi + S2 + 2s3 - 4M2) - l(si + S2- 553) + 
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■]:Ml[J{s^) + J{h)] - ^{si - S2)[J{h) - J{h)] + \sj{h) 



1 



967r2 



2^1 (S3 



2M2) - + S2 - S3 - 4M2) + 6r4M2 



1 



+h{.si + S2- 2M^) + -(si + S2 - 7s 



10 



Mi 



4(^1 - S2)[J{h) + J(S2)] - ^(Sl + S2 - 5M2)[J(si) - J{h)] 

o o 



1 



967r2 



-2/2(51 - S2) + -(si - S2) 



(44) 



where Sj = Si/M"^^ and the Tj are hsted in Tab. 1. 

As for A^TT its expansion is now known up to the two loop level ||3^. Here we 
need only the first two terms of the expansion which are known since a long time 
13^,111: 



s-Mi 



FlAfl{s,t,u) 



1 

12 



[{t - u)^ - 2Mls + 4M^] [ J(t) + J{u)] 



~^{t - u){s + 2Ml)[J{i) - J{u)] + \{s' 



Mt)J{s) 



u? + s^] 



-3hM^ + l2hMl{s - Ml) - -{t - n) 



6 



(45) 



The form factors F and G as defined in ( ^OD can be easily reconstructed from A^^^^ 
and the corresponding barred functions at each given order, via the simple relation: 



R 



pn ^(n) 



TT TTTT 



Ml - Q2 



R = F,G . 



(46) 



One could wonder whether it is possible to experimentally disentangle the contri- 
bution of A^TT to some of these form factors. We are convinced that this is not the 
case. The reason being that A^^^ contributes only to Fs, which is very difficult to 
measure, and that there are in addition other contributions to Fs- 
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Table 2: Integrated branching ratios for yQ'^ < Qmax predicted with CHPT at a 
given order 0{p^) and from vector meson dominance models (VMD), see text. 



mode Qma.[MeV] 0{p^) 0{p^) 0{p^) VMD 



400 4.22 X 10"^ 6.66 x 10"^ 7.34 x 10"^ 7.40 x 10"^ 
2'KUr 500 1.57 X 10"^ 2.92 x 10"^ 3.45 x 10"^ 3.72 x 10"^ 
600 3.40 X 10"^ 7.57 x 10"^ 9.73 x 10"^ 1.20 x 10"^ 



500 2.00 X 10"^ 4.19 x 10"^ 4.19 x 10"^ 

37rz/^ 600 4.31 x 10"^ 1.07 x 10"^ 1.23 x 10"^ 

700 2.21 X 10"^ 6.55 x 10"^ 9.51 x 10"^ 



5 Numerical Results for the Branching Ratios and 
the Structure Functions 

5.1 Two Pion Decay 

At first, let us discuss integrated branching ratios ^^2-K{Qmax) 
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BR2.(g„a.) = BR, /"™"^ rfQ2,i^(g2) (47) 

The results from CHPT are given in Tab. ^. From the convergence of the expansion 
in p^ we conclude that the CHPT expansion truncated at this order works fine up 
to Q^ax = 500 MeV. 

We also compare with the prediction from a VMD model, viz. from the model 1 of 
1^, which has been implemented in TAUOLA |Q. This VMD model parameterizes 



Fy in terms of a coherent superposition of a p and a p' Breit-Wigner, with an 
overall normalization fixed by matching to the O(p^) chiral prediction. It gives a 
good parameterization of e~^e~ — * 27r annihilations in the range covered by the tau 
mass. We find that in the range up to 500 or 600 MeV, where we trust CHPT, the 
predictions of CHPT and from the VMD model agree well. 
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Figure 1: Differential decay rate for r 2-ni'T-: Predictions by CHPT at Oiji^) 
(dashed), at 0{p^) (dashed-dotted), at 0{p^) (solid), and from a vector meson dom- 
inance model (dotted), compared with experimental data from CLEO (dots with 
error bars) 




[GeV^] 



In Fig. |1] we plot the two pion invariant mass spectrum normalized to the elec- 
tronic branching ratio of the tau. We plot the predictions from CHPT, together 
with the prediction from the VMD model in [H] and preliminary data from CLEO 



3J] . Note that we used a simplified approach to fix the overall normalization of the 



data in [Q. We multiplied the spectrum from with a normalization factor N 
and determined by fitting the data to the VMD prediction of . Of course, the 
normalization should instead be taken from the data. In fact we suggest a careful 
reanalysis of the low energy part of the spectrum and its absolute normalization in 
order to compare it with the CHPT prediction. 

It is of some interest to understand how sensitive r 2txVt- decays are to the 
pion charge radius {r'^)vi which is defined from the expansion of Fyi^Q"^) in terms of 

FviQ') = 1 + lir'rvQ' + clQ' + (g^) + Om (48) 
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where fy{Q'^) is given in ( |36D and is very small numerically. According to the pre- 
vious results, we consider this expansion valid up to Qmax = 500 MeV. Furthermore 
we neglect theoretical uncertainties due to higher order corrections and assume that 
Cy is known exactly. For a discussion of these points see Sec.^ . 

Given a number N of events r — >■ 2711'^- with hadronic invariant mass squared 



in the interval 4M^ 



where 



with 



Qmax^ the precision with which (r )y can be measured is (see 




df{x;py 
dp ^ 



dx 



1 



-1/2 



N 



7.37 fm^ 



(49) 



(50) 



piQ 



2, 1 dT^AQ^) 



T, dQ^ 



R 



4Af2 



dQ^p{Q^ 



So given A^^ decaying taus and an detection efficiency r], we have 
N = N^T] BR(r ^ 2niy^, < 0.25 GeV^) 



(51) 



(52) 



where the branching ratio BR(r — > 2711/^-, Q < 0.25 GeV ) = 3.6 x 10 according 
to Tab. I 

Based on this, we now give rough order of magnitude estimates for the possible 
statistical accuracy of present and future experiments. We assume an efficiency of 
r] = 30%. CESR has at present about 5 x 10^ taus, thus the possible statistical 
accuracy is of the order of 0"(r2)^ = 0.1 fm^. A b-factory might have 5 x 10^ taus per 
year. Assuming 3 years of running time, this leads with the assumptions mentioned 
above to a possible statistical accuracy of the order of cTij^z)^ = 0.02 fm^. 

These numbers have to be compared to the accuracy evaluated from vre scatter- 
ing. With the assumptions we made above, the present result is {r'^)y = 0.431 ± 



0.010 fm^ [0]. So it seems difficult for tau decays to become competitive with ire 
scattering for the determination of {r'^)y- Nevertheless, tau decays can provide an 
interesting cross-check. 



5.2 Three Pion Decay 



The numerical results for the integrated branching ratio with < Ql^iax CHPT 
at 0{p^) and O(p^) are given in Tab. 0. The series expansion of CHPT does not seem 
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Figure 2: Integrated structure function Wa{Q'^) for 2-n~Ti'^: CHPT prediction at 
0{p^) (dashed), 0{p^) (solid) and from a vector meson dominance model (dotted). 
The four functions wa and Wc for both modes 2i:~'n^ and 27r°'7r~ all look identical 
within the resolution of this diagram. 



0.6 — 

0.4 — 

0.2 — 

Eq 




0.2 0.25 0.3 0.35 0.4 

[GeV^] 



to behave very well. Even for Qmax = 500 MeV, the integrated one loop contribution 
is already around 45 % of the tree level result. This means that to have a CHPT 
prediction with a reasonably small theoretical uncertainty for such a quantity one 
should stop at Qmax Well below 500 MeV, and this would be very difficult to test 
experimentally, because of the strong phase space suppression of this region. On 
the other hand we observe that there is a fair agreement with the VMD model 
numbers up to Qmax = 600 MeV, which is better than what happens in the two pion 
case. This can be understood in terms of the fact that in the two pion case, the 
p(770) resonance is very close, whereas in the three pion case, the nearest three pion 
resonance, the ai(1260), is much further away. Whether the VMD model is a good 
representation of the experimental data even at such a low energy though, has still 
to be verified. 

Beyond the integrated decay rates, this decay mode has a very rich structure, 
which in principle could be investigated in detail experimentally. In fact in [27 it has 
been shown that it is possible to extract all three form factors from a measurement of 
all angular distributions. One could then compare the measured form factors to the 



21 



Figure 3: Integrated structure function Wc^Q"^). CHPT at O^p'^) for 2n (solid) 
and for (dashed-dotted) and form the VMD model (dotted, identical predic- 

tion for both charge modes) 

^^[GeV^] 
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analytic formulae of CHPT. This extraction is however very difficult in practice, and 
moreover the phase space suppresses considerably the region where the comparison 
makes sense. For these reasons we find more useful to show directly the curves for 
the integrated structure functions wx in Figs. |^-|^. We have chosen six out of the 
nine structure functions which are present, mainly because the missing three are too 
difficult to be measured and not particularly interesting. Note that, in accordance 
with our definitions in Eqs. (^,^), the Cabibbo angle cos^ 6'c is factored out from 
hadronic matrix element, and so the structure functions wx do not include this 
factor. 

In Fig. we plot WAiQ'^) for the 27r~7r^ final state. It turns out that Wyi(Q^) 
and wciQ"^) for both charge states, 2n^n^ and 2tc^7c~ are all very similar, so the 
four plots can not be distinguished from each other. The equality of the structure 
functions for the two charge modes suggests that wa and wc are dominated by 
the [210] partition. We have explicitly verified that this is the case: for example 
at = IhMl the [210] partition contributes 6.73 x 10"^GeV"^ to wa for both 
modes, whereas the [300] partition contributes 2.21 x lO^^GeV"'' for the 27r~7r'^ and, 
according to Eq. (^), one quarter of this for the 27r°7r~ mode. The interference of 
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Figure 4: Integrated structure function We{Q'^)- CHPT at 0{p^) for 27r tt^ (solid) 
and for (dashed-dotted) and from the VMD model (dotted, identical predic- 

tion for both charge modes) 
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the two partitions is completely negligible. 

The reason why wa ~ wc at low energy is that 

Wc = Wa - 2xl\Fi - F2\^ 

= WA-8xl\Hf. (53) 

Since H is zero at tree level, the difference vanishes at leading order. Moreover 
this difference starts as the square of a quantity of OijP'). Using the language of 
the sixties we can say that the vanishing of this difference near threshold is a Low 
Energy Theorem (LET), and that it receives corrections only at next-to-next-to- 
leading order. We also notice that if is a function antisymmetric under exchange 
of Si and S2, so that its modulus squared has a zero along the line si = S2- At low 
energy, where the distance between this line and the boundaries of integration is not 
large (in units of M^), the presence of the zero produces an additional suppression 
of the integral. The similarity of these two functions near threshold was already 
found in Ref. [^, in the framework of a VMD model. Here however, we can 
give a detailed algebraic account of why this happens. For all these arguments we 
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Figure 5: Integrated structure function WsAiQ"^) for 27r tt"*". CHPT prediction at 
0{p^) (dashed) and 0{p^) (solid), and from a VMD model (dotted). 
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believe that this fact should be verified experimentally even at energies well above 
those where one would trust a one loop CHPT calculation. 

In Fig. H and § we plot WDiQ"^) and we{Q'^)- These are much smaller than wa 
near threshold, and certainly more difficult to be measured. However they seem 
rather interesting from a theoretical point of view. First of all for both of them 
the 0{p^) contribution vanishes. Secondly, the O(p^) predictions from CHPT differ 
strongly for the two charge modes. According to the discussion in Sec. |3.3.2| , this fact 
suggests that near threshold these two structure functions are strongly influenced by 
the partition [300]. Using the decomposition of the form factors given in Eq. (pSj), 
we have verified that this is exactly what happens. As can be seen in Tab. I the 
change of sign in the interference contribution is responsible for the change of sign 
of the whole integrated structure functions in the two different charge modes. As 
far as we know up to now there are no data for these structure functions so close 
to threshold: so this is a real prediction of CHPT. Moreover this change of sign is 
absent in all the models of which we are aware which have been used to describe 
this decay channel. It is then important to ask how reliable this prediction is and up 
to what energy it can be trusted. These questions are especially difficult to answer 
here because we have only a leading order calculation (since the O(p^) contribution 
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Figure 6: Integrated structure function Wsa{Q'^) for 27r°7r . CHPT prediction at 
0{p^) (dashed) and 0{p*) (solid), and from a VMD model (dotted). 
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[300] 


[210] 


interf. 


total 


(27r-7r+) 
(27r°7r-) 


8.63 X 10~^ 
2.16 X 10-^ 


1.02 
1.02 


-1.56 
0.779 


-0.450 
1.82 


(27r-7r+) 
(27r°7r-) 


4.62 X 10-^ 
1.16 X 10-^ 


0.206 
0.206 


2.08 
-1.04 


2.29 
-0.837 



Table 3: Splitting of the integrated structure functions we and wd into the contri- 
butions of the partitions [300], [210] and their interference, at = 15 M^. 
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Figure 7: Integrated structure functions ly^^lQ^)- Prediction from CHPT for 27r tt 
(solid), 2'7r°7r~ (dashed-dotted), and from a VMD model for 27r~7r^ (dotted), 
(dashed) 
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vanishes). The only thing we can do is to check how sensitive the prediction is to 
the values we use for the low energy constants. We concentrate here on the [300] 
state, since it is the one responsible for the new effect: 



[300] 



iV2 



727r2 

4 , 



- si + Ml 



+ -J[s,){Mi-s,) 



2 _ 

g^(S2)(S2 



Ml 



2 , 



9 



J{s,){s,-Ml). (54) 



First of all we notice that this partition has no contribution from tree level. This is 
not the case for the [210] partition, which starts as: 



i\/2 



[210] 



3+0(p^) 



(55) 



This explains why the interference contribution is much bigger than the one from 
the [300] partition alone. Then we may easily see from the definition of these two 
structure functions given in App. [A.3| that the W£, is mainly sensitive to the real 
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Figure 8: Dalitz plot distribution of the 2-k -k^ final state in si, $2 for = 
0.36 GeVl 

Si [GeV^] 
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part of F^'^'^^ , whilst We to the imaginary part. This means that only the numbers 
for wd may depend on the combination of low energy constants which occurs in 
^|300] rjj^g value used for the combination li + I2 — 3 in the numerical calculations 
is 1.4, resulting from the central values given in Tab. |I|. The uncertainty on this 
number can be estimated from the same Tab. |I| to be ±1.1. So the possibility that 
the real value of this combination be twice as much, or very close to zero cannot be 
excluded. We have checked that by changing the value of Ti ± [2 — 3 by ±1.4 changes 
wd in the 27i~7i~^ mode by ±4.4 x 10"'', at = 15M^. As expected we remains 
practically unchanged. Even in the worst case for wd, however, there remains a 
sizable difference between the two charge modes. The effect of higher orders in the 
chiral expansion and the related questions of how far in energy one can trust this 
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prediction will remain unanswered until one calculates the form factors at two loops, 
which is beyond the scope of our analysis. We may just argue that since there are 
no low-lying resonances contributing to this particular three pion state, one of the 
possible sources of large higher order corrections is excluded. 

All in all we can say that CHPT predicts a sign difference in the two structure 
functions W£, and we for the two charge modes near threshold. This prediction will 
have interesting consequences in the next subsection, where we compare with VMD 
models and available data. 

In Fig. ||we plot WsAiQ"^) for the 2'k~'k^ state. The corresponding curves for the 
27r°7r~ state in Fig. ^have a very similar shape but a different overall normalization. 
In fact, the ratio w^'^ I'^'s'a^ is close to 4, which according to Sec. |3.3.2| implies 
that the scalar form factor at low energies is dominated by the [300] state. 

We find that the structure function wsa is very small compared to wa, so the 
spin-0 contribution to the integrated rate is negligible. The only chance to measure 
the spin-0 part is via its interference with the spin-1 part in the structure function 
wsb{Q'^), which is plotted in Fig. ^ 

In Fig. § we consider the Dalitz plot distribution dT / (dQ^ dsi ds2) in si, S2 for 
fixed = 0.36 GeV^. We display the O(p^) prediction for the 27r^n^ mode. At 
0{p'^), the Dalitz plot density depends only on si + S2. As seen from the figure, 
this feature seems to persist at O(p^). Let us analyze this issue quantitatively. To 
describe the Dalitz plot, we define s+ and s_ by 



si - S2 



(56) 



and then replace s+,s_ by dimensionless variables x,y with < x < 1 and —1 < 
y <1 via 



(57) 



where 



- ml 



\ 



[s+{s+ - Ami) ~ 4m2((52 _ ^2)][g2 - s+ - ml] 
+ 3m2 - s+ 



(5J 



28 



We perform a least-square fit to the Dalitz plot density as predicted by CHPT, using 
a fit function 

Pq2 (x, y) := , — = a{x + b + cx^ + dx^ + ex"^ + fx{x - 1)?/^) (59) 

a(y^ dsi ds2 

Note that this is a reasonable ansatz for the y dependence. Firstly, we must have 
Pq2{x, y) = Pq^^x, —y) because of Bose symmetry. Secondly, at x — and at x — 1, 
^max ^ g gQ dependence on y"^ must go to zero at x = 0, 1. 

We choose to discuss the CHPT predictions for = 0.36 GeV^. At this value 
of Q^, the result from the fit is 

a = (3.0360 ± 0.0095) X 10"^^ GeV"^ 

b = (5.251 ±0.039) x 10"^ 

c = (0.786 ± 0.023) 

d = (0.190 ±0.053) 

e = (0.189 ±0.036) 

/ = (0.0645 ± 0.0026) (60) 

This fit never deviates from the CHPT prediction by more than 2%, with an average 
deviation of less than 1%. It is seen that the y'^ dependence is very small. In fact, 
taking into account that the coefficient x(x — 1), which multiplies y^, is less or equal 
to 1/4, we find that the leading y dependence does not exceed 2%. We have checked 
that, as one would expect, the y'^ dependence is even smaller for smaller Q^. 

This Dahtz plot distribution, as predicted by CHPT for < 0.36 GeV^, differs 
strongly from the behavior at high in the resonance regime, where the p reso- 
nances lead to pronounced structures in si and S2 (resonance bands for fixed si or 
^2) 



5.3 Comparison with Vector Meson Dominance Models 

In this subsection we will compare with the low energy behavior of the phenomeno- 
logical models in 1^, 1^, §]. The simplest VMD model which one can build for this 
channel (see Ref. 0) is based on the decay chain r — aiz/^, ai pn, p — >• 27r and 
a transverse ai propagator. In this case the amplitude contains only a spin 1 part 
{Fs = 0), and the three pions are only in a [210] partition state (Ff°°^ = 0). The 



comparison to the data ||3^, ^ shows that the model works well, which means that 
the assumptions made are reasonable. 

However it is clear that these assumptions need not be strictly true in the physical 
reality, so the authors of Refs. [0, || have tried to include in the VMD model a 
nonzero scalar form factor. Two possible sources for a non-vanishing Fs, are a 
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pseudoscalar three-pion resonance, the vr', or the non-transverse component of the 
off-shell ai propagator. A model for the vr' contribution is given in . The numerical 
impact of the vr' depends on a parameter f.,^/. In [0 and in its implementation in 
TAUOLA a particularly large value from 0] has been chosen, which is probably 



several orders of magnitude too large |^ . However, we find that even with this high 
value of /tt', the contribution from the vr' to the scalar form factor at low energies 
is much smaller than the predictions from CHPT. This indicates that — at least at 
small energies — there are additional contributions to Fg. The off-shell contribution 
of the ai to Fs is discussed in A specific model is constructed by matching to 
the 0{p^) prediction of CHPT (including the ^ effects), and we will compare 
this model with the CHPT predictions. 

On the other hand, the presence of a [300] component in the spin 1 form factor 
has never been proposed in any of these models. With our calculation we can make 
a detailed analytical comparison between VMD models and the CHPT amplitude 
at low energy. The main conclusion is that even in the region close to threshold the 
VMD model works rather well. 

First we consider the structure functions which only involve the Fi. After a 
proper normalization of the form factors in the VMD model, which takes into account 
the CHPT expressions at O(p^), the agreement for the spin one spectral function 
at low energy looks very good. This can be seen in Fig. |^, where we plot wa, i-e. 
the main contribution to the total decay rate. Moreover the [300] component of Fi 
starts only at one loop, which means that in the chiral expansion this is algebraically 
suppressed with respect to the [210] part. 

However, the - numerically rather small - structure functions Wo and We are 
sensitive to the [300] part via interference with the [210] part, as we have seen in 
the previous subsection. Comparing the VMD model with the CHPT prediction for 
these structure functions in Figs. |^ and we find good agreement for 27r°7r^, but 
flat disagreement for 27r~7r^. However, at larger Q^, experimental data for wd and 
We in 27r^7r^ are available, which agree with the VMD prediction [^8[ Note 
that the left-right asymmetry A^ji^Q'^) measured by ARGUS is proportional 
to we{Q'^)/wa{Q'^) [0, and confirms the VMD prediction in sign and magnitude 
down to = 0.8 GeV^. We conclude that, unless the higher orders in the chiral 
expansion completely change the CHPT result, somewhere between threshold and 

= O.SGeV^ there must be a zero for both structure functions in the 27i vr^ mode. 
It would be extremely interesting to verify this zero experimentally, or, as a minimal 
option, to verify the existence of a difference between the two charge modes. 

Next we consider the scalar form factor: Although Fs is nonzero already at tree 
level, it is kinematically suppressed: 

^ mho,, 
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Though Q^,S3 and are all counted as quantities of order p^, so that the ratios 
M^/Q'^ and s^/Q"^ are algebraically of order one, it is clear that numerically they are 
smaller than one. For example at threshold Fs = iV2/F^{— 1/24) = —1/16 x Fi. 

In Figs. ^,11 and 0, where we plot the two structure functions wsa and wsb, one 
can see the comparison of the VMD model in § to CHPT at 0{p^). We find that 
the VMD model, which by construction reproduces the 0{p^) prediction, does not 
reproduce well the O(p^) prediction from CHPT. The reason for this discrepancy 
(which is larger for wsa, since it contains the scalar form factor squared) lies in the 
fact that in the model in the [300] part of Fs is added as a constant, without a 
resonance factor enhancement. 

Summarizing our comparison of CHPT and VMD predictions at low energies, 
we have found that VMD gives a good description of the dominant spin-1, [210] 
contribution. However, CHPT shows that both Fs and fI^^^\ though small, are 
not exactly zero. This fact, which is a prediction of CHPT, has still to be verified 
experimentally. Our analysis shows that the best place to look for these parts of the 
amplitude, is where they can interfere with the "big" components: the presence of 
Fs should be detected by measuring a nonzero wsb, whereas the presence of 
should be discovered by measuring a sizable difference (possibly a sign difference) 
between the two charge modes for wd and we- 



6 Summary and Conclusions 

Chiral perturbation theory (CHPT) provides model-independent predictions for 
hadronic matrix elements in the low energy region below 500 ■ ■ ■ 600 MeV. It does 
not contain additional assumptions beyond the fact that the strong interactions are 
described by the QCD Lagrangian and that QCD possesses an approximate chiral 
symmetry which is spontaneously broken. We have evaluated tau decays into two 
pions (and tau-neutrino) to two loops and decays into three pions to the one loop 
level. The branching ratio into the phase space region with small enough invariant 
hadronic mass for CHPT to be applicable was found to be about 4 x 10~^ for the 
two pion mode and about 10~^ for the three pion mode. And so the predictions of 
CHPT for r — > 27rz/^ are testable at present machines (LEP and CESR), while in 
the case of r ^ SttUt-, future facilities with very high r production rates seem to be 
required (b factories, r-charm factory). 

In the case of the decay r — >■ 27?//,-, the predictions for the invariant mass dis- 
tributions can be tested at present experiments, and the spectrum can be used to 
extract the pion charge radius (r^)^. 
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As for the decay t Svrz/^, a detailed comparison to the CHPT predictions 
for the form factors near threshold requires very high statistics, mainly because of 
the phase space suppression. For this reason we have tried to identify a few spots 
where the consequences of the approximate chiral symmetry of QCD can be tested 
experimentally with a reasonable statistics. These are: 

1. the very close similarity of wa and wq near threshold, which seems to extend 
well beyond the very low energy region; 

2. a zero and change of sign in wd and We between threshold and ~ O.SGeV^ 
in the channel 27r""7r''", or at least a difference between the two charge modes 
near threshold. This would be the first evidence for the presence of the [300] 
partition state in this decay channel; 

3. the presence of a scalar form factor as predicted by CHPT, to be detected by 
measuring wsb- 

Comparing our results to predictions from vector meson dominance models, we 
find overall a reasonable agreement in the low energy region. The structure function 
Wa, which dominates the decay rate, is described well by VMD models. We have 
found, however, some interesting discrepancies in certain (numerically rather small) 
structure functions. These discrepancies are related to the [300] partition state and 
to the scalar form factor, which both are missing or underestimated in vector meson 
dominance models. 
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A Appendix 

A.l General conventions 

Consider the decay of a tau into n pions, 

r(p, S) Vriq, s')'Ki{pi)'K2{p2) ' " " T^n{Pn) (62) 

Here s and s' denote the polarization four-vectors of the tau and the neutrino, 
respectively. We define the total hadronic momentum Q by 

Q = Pi + P2 H h Pn (63) 

and in the case of three pions, we will use Dalitz plot invariants si, S2 and S3 defined 
by 

S1 = (P2+P3)' (64) 

and cyclic permutations. 

The differential decay rate dVn is given by 

dTr. = ^\M\^d^n (65) 

where the phase space element d^n is 

dt„^(2.mQ+,-p)p||^np|l|g^ (66) 

and 

M = ^cosOcL^H^ 

= {7ri{pi)n2ip2)...7iniPn)KiO)-A^{0)\0) (67) 

A. 2 Decomposition of the form factors in terms of partition 
states 



As we discussed in Sect. p.3.2| , we have three possible partition states, two belonging 



to the class [210] (which we will indicate as [210]a and [210];,) and one belonging 
to the class [300]. The matrix elements of these states are easily constructed by 
performing the appropriate symmetrizations and antisymmetrizations with respect 
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to momenta exchanges, as described in Sect. p.3.2| . The result can be expressed in 
terms of the three form factors F, G and H, and reads as follows: 



^[210] a 
Q[2W]a 
^[210] a 



^[210]6 
Q[2l0]t 



G23 + H- 



23 



G12 — - (-^23 + ^23 



H12 — - {F23 



H23) 
G23 — H23] 



1 



G 



23 



77 [—-^13 + -^23 
O 

7 [-^13 + -^23 ~ ^"13 
O 



Hi3 — i/23] 

G\3 + G23 ~ -f^l3 + -^^23] 
G23 ~ -f^l3 ~" -^^23] 5 



(6J 



(69) 



^[300] 
(^[300] 
^[300] 



[F12 + Gi3 + G23 — — H- 



23j 



G12 + - {Fi3 + F23 + Gi3 + 6*23 



-f^l3 + 



23 y 



H12 + - ( — -^13 + -^23 + ^13 ~ ^23 + — H23) 



(70) 



where Xij = X{si, sj, Sk), X = F,G,H and k ^ i ^ j. To reconstruct the form 



factors in the two charge modes one has to use the relations: 

X12 = x[^°°i + X'^^^l" + X'^^^l" 



X 



(+) 
12 



2J^[300] _ j^[210]„ _ j^[210]i, 



(71) 



A. 3 Structure functions for the three pion final state 



In this section we will briefly review the formalism of structure functions for the 
decay of the tau into three pions, and display formulae relevant for the present 
paper. For more details, the reader is referred to [^, ^ . 



The decays are most easily analyzed in the hadronic rest frame 

= g=pi+pl+pl (72) 

The orientation of the hadronic system is characterized by three Euler angles a, f3, 
7, as introduced in They can be defined by (0 < a, 7 < 27r; < /3 < tt) 



COS a 



[riL X n_L) 



Fl X nrWriL X nj_| 



sm a 



Fx, X UrWriL X nj_| 
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COS 7 



sm7 



fiL X n± 



cos P = Hl ■ fi± 



(73) 



where ul is the direction of the laboratory in the hadronic rest frame, n± = {pi x 
P2) /(|pi X P2I) is the normal to the pion plane (here we assign the momenta according 
to 1^2! > \pi\), nr is the direction of flight of the r in the hadronic rest frame and 
ns = Ps/IpsI- 

From these definition it is obvious that /3 and 7 are observable even if the tau 
rest frame cannot be reconstructed, whereas a does require this knowledge, a could 
be measured at a tau-charm factory where the tau pairs are produced almost at rest 
and therefore the tau rest frame is known, or if the tau direction can be measured 



with the help of vertex detectors [p]] . 

The Euler angles also serve to parameterize the phase space 



1 1 - Q2 (Iq2 dad-id cos (3d cos 9 

-dsi ds2- 



(27r)5 64 M2 Q2 



27r 27T 



(74) 



where 6 is the angle between the direction of flight of the r in the laboratory frame 
and the direction of the pions as seen in the r rest frame, and cos^ can be calculated 
from the energy of the pion system with respect to to the laboratory frame and 
beam energy Ebeam 

2E,/E,M'^ - M2 - Q2 

cos 6' = , (75) 

(M2-Q2)^l_M2/i?2_ 

The contraction of the leptonic and hadronic tensors can be expanded in a sum 



(76) 



X 



In general, H^'^ can be characterized by 16 independent real functions. In our case of 
a three pion final states, there are restrictions due to G parity and Bose symmetry, 
which leave 9 independent functions. In a convenient basis they are given by 



Wa = {xj + xl) + {xj + xl) iFsp + 2{xiX2 - x^) Re (FiF^) 

Wb = xl\Fs\^ 

Wc = ixj-xl)\F^\^ + {xl-xl)\F2\^ + 2{xiX2 + xl)Re{FiF;] 

Wd = 2[xiX3\Fi\'^ - X2X3\F2\^ + X3{x2- xi)Re{FiF;] 

We = -2x3{xi + X2)lm{FiF;) 
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Wf = 


2x4 [xi Im (F1F3*) + X2 Im (F2F3*)] 


Wg = 


-2x4 [xi Re (F1F3*) + X2 Re (F2F3*)]] 


Wh = 






Wi = 


-2x-^xa \ Re{F^Fn ) - 


Re fFoFo*)! 


yVsA — 






WsB = 


2^02 [xi Re (FiF;) 


+ X2 Re (F2F;)] 


Wsc = 


-2^[xiIm(FiF4*) 


+ X2 Im (F2F4*)] 


WsD = 


2^02x3 [ Re (FiF;) 


- Re(F2F;)] 


WsE = 


-2^X3[Im(FiF;) 


- Im(F2F;)] 


WsF = 


-2^^Xi \^[F^Fl) 




WsG = 


-2Jq^x^Re{F^Fl) 





(77) 



The variables Xj are defined by 



Xi 

X2 
X3 

X4 



<fl-<lt 

0.2 - O3 

V 

Oi 



V 

-02 



iQ'^x^ql 



(7J 



where (gf ) denotes the x {y) component of the momentum of meson i in the 
hadronic rest frame. They can easily be expressed in terms of si, S2 and S3 |27 |. 

The hadronic structure functions Wx depend only on Si, S2 and Q^. The corre- 
sponding leptonic Lx depend on the Euler angles a, /3, 7 and on E^. The relevant 



formulae can be found in 27 



The structure functions can be measured by observing angular distributions in 
the Euler angles /?, 7 and, if the tau rest frame is known, a, and taking moments 
(m) with respect to products of trigonometric functions 



(m) := 



2(M2 - g2 



d cos P dj 
2 2^ 



(79) 



As shown in measuring suitable moments involving (3 and 7 only (and, in 
fact, an energy ordering sign(si — S2) in some cases to avoid vanishing due to Bose 
symmetry) allows to extract all the individual structure functions except for Wsc 
and WsE- If the tau rest frame and hence a is known additionally, Wsc and Wse 
can be measured, too. 
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